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ABSTRACT 

We have discovered a correlation between the observed peak spectral energy 
-Epk.obs and the Euclidean value of < V/V ma _ x > of gamma-ray bursts (GRBs). 
We present the evidence for the correlation in the GUSBAD catalog and use 
it to derive the luminosity function of GRBs without using any redshifts. The 
procedure involves dividing GUSBAD GRBs into five spectral classes based on 
their -E p k,obs- The overall luminosity function is derived assuming that each of 
the spectral classes contributes a gaussian luminosity function. Their central 
luminosity is derived from the observed Euclidean < V/V max >. We explore 
various forms for the GRB rate function GR(z) in predicting redshift distributions 
of GRBs detected by Swift. We find that GR(z) peaks at a higher redshift 
than the typical star formation history currently favored in the literature. We 
consider two examples of GR(z) that sucessfully predict the observed redshift 
distribution of Swift GRBs. With the luminosity functions in hand, we convert 
the £7pk,obs ~ V/Vjbbx correlation into an -Epk,obs — -^iso correlation and a rest frame 
Epk — L iso correlation. In comparing the — L mo correlation with a published 
correlation based on GRBs with known i? p k,obs and redshifts, we discuss the effect 
of Malmquist bias. 

Subject headings: gamma rays: bursts 

1. Introduction 

The GRB luminosity function plays a central role in the interpretation of the observed 
source counts and redshift distributions of GRBs detected with different instruments in dif- 
ferent energy bands and to different detection limits. Of particular interest are the resulting 
rates and the variation with redshift, since these must reflect the properties of the GRB 
progenitors. 
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Ideally, the luminosity function is derived from a well defined flux-limited sample of 
GRBs with redshifts. For most GRB surveys, it is not possible to collect such a sample with 
well defined gamma-ray and optical flux limits. We will discuss the situation for the Swift 
mission in Section 4. 

In the absence of usable redshift sta tistics, most studies of the lumin osity function use 
l uminosity criteria, such as var iability (jFenimore fc Ramirez-Ruizl 120001 ) and spectral lag 
(INorris. Marani fc Bonnellll2000l ). We have reviewed several luminosity functions resulting 
from these studies and found that they generally are not compa tible with the s ource counts 
and the value of < V/V max > for GRBs in the BATSE catalog JSchmidtl 120041 ) . As long as 
this is the case, they cannot provide reliable predictions about redshift distributions, etc. 

The correlation between the rest fr ame spectral peak energy E p ^ and the isotropic- 
equivalent radiated energy E iso found by lAmati et al. J2002h seemed to be very promising 
as a l uminosity indicator. However, it has been shown (INakar fc Piranll2005l ; iBand fc Preece 
20051 ) that for a large fraction of BATSE b ursts there is no redshift that satisfie s the 
Amati relation ; for fu rther discussion al so see iGhirlanda. Ghisellini fc Firmanil (120051 ) and 
Nakar fc Piranl (120051 ). According to|Lj (120071 ). a redshift degeneracy in the Amati relation 
makes it impossible to derive redshifts larger than 0.9 with usable a c curac y. Using durations 
and spectral parameters for GRBs detected by Swift iButler et al.l (120071 ) find a E-^ — Ei so 
correlation that i s inconsistent with the Amati relation at > 5cr significance, with double 
the scatter. Both lLloyd et al.l (120001 ) and lButler et al.l (120071 ) discussed the effect of detector 
thresholds on the derivation of E p ^ — E iso correlations. The latter authors argue that all 
pre- Swift E v ^ — E [so correlations are likely unrelated to the physical properties of GRBs. 



In a study of the GRB luminosity function from 82 HETE-2 GRBs (IPelangeon et al. 



20081 ). an Amati- type correlation was used to derive pseudo-redshifts for the 62 GRBs lacking 
an observed redshift. This study provides information about low-luminosity GRBs and X-ray 
flashes that are not represented in the BATSE based GUSBAD catalog. 

In the present study, we introduce a new correlation between the peak spectral energy of 
the observed peak flux i^obs an d the Euclidean value of < V/V ma _ x >. The correlation was 
discovered for GRBs in the GUSBAD catalog 0. Since the Euc lidean value of < V/V max > 
is a geometric cosmological distance indicator (jSchmidtl boOll ). the peak luminosity of a 
subgroup of GRBs with given £7 p k,obs can be derived without knowing any redshifts. Using 
the -Epk,obs — V/V max correlation is not subject to any of the problems mentioned above for 
the Amati relation. 



1 Available at http//www. astro. caltech.edu/~mxs/grb/GUSBAD. 
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The approach in our paper is the following. We use GUSBAD fluxes in the four BATSE 
DISCLA channels to divide them into five spectral classes based on their i? p k,obs values. We 
assume that the luminosity function of each spectral group is a gaussian of log L with a value 
of <7\ gL large enough so that the overall sum of the five luminosity functions is reasonably 
smooth. In addition, we assume that the overall rate density of GRBs varies with redshift 
as GR(z). The central luminosity logL c of each spectral gaussian is derived by a process 
of iteration, until the corresponding value of V/V max for the spectral class agrees with that 
observed. 

The procedure allows deriving the overall GRB luminosity function from the GUSBAD 
catalog without using any redshifts. The main free parameter is the rate function GR(z). To 
evaluate the predicted luminosity and redshift distributions for various forms of GR(z), we 
do need a sample with redshifts. We evaluate and use Swift detected GRBs with observed 
redshifts. For the successful luminosity functions, the process yields a rest frame E pk value 
for each of the spectral classes. This leads to -E p k,obs — -^iso and E pk — L iso correlations entirely 
based on data in this paper. 

The use of V/V max and the derivation of £7 p k,obs from GRBs in the GUSBAD catalog 
are discussed in Section 2, as well as the evidence for the -E p k,obs — V/V mSuX correlation. 
The framework for deriving the luminosity function is discussed in Section 3. The status 
and completeness of Swift observations is covered in Section 4. In Section 5 we present 
two luminosity functions with differing redshift dependence GR(z) and test them on Swift 
redshift and luminosity distributions. The £pk,obs ~~ -^iso and E pk — L iso correlations are 
derived in Section 6. A discussion and summary follows in Section 7. 



2. The E pk}0hs - V/V max Correlation for GUSBAD GRBs 
2.1. Euclidean Values of V/V max 



The GUSBAD catalog flSchmidtl 120061 ) is based on observations with the BATSE LAD 
detectors which provide output in four energy channels, viz. 20 — 50 keV (ch 1), 50 — 100 
kev (ch 2), 100-300 keV (ch 3), and > 300 keV (ch 4). The catalog lists peak p hoton fluxes 
for ch annels 2 and 3 together. These were derived assuming a Band spectrum (IBand et al. 
19931 ) with a = —1.0, (3 = —2.0, and E = 200 keV. For the present study, we have derived 
for each GRB peak photon fluxes in each channel based on the two brightest illuminated 
LAD detectors. 



At the outset of this study, it was not clear whether the four channels of the BATSE 
LAD detectors could be used for low-resolution spectrophotometry. We decided to use only 
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GRBs with a 50 — 300 keV peak photon flux P exceeding 0.50 ph cm s , i.e. twice the 
effective photon flux limit of the GUSBAD catalog. This still leaves a large sample of 1319 
GRBs. The Euclidean value of V/V ma _ x is simply V/V maiX = (P/Pu m ) 
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The value of < V/V max > for uniformly distributed objects having any luminosity func- 
tion in Euclidean space is 0.5. This value derives from the fact that volumes are ~ i? 3 and 
areas ~ R 2 . This will not apply in non- Euclidean space, so for cosmological objects the 
Euclidean value of < V/V ms ^ > will deviate from 0.5, the more so the larger the typical 
redshift, i.e., it is a cosmological distance indicator. 

Concern has been expre ssed in the literature about the use of VjV ma x particularly 
when sample thresholds vary (jBandlll992l ; |Petrosianlll993l ; lHartmannlll993l ) . The problem of 
varying threshold can be handled if each object has its own value of Pu m as in the GUSBAD 
catalog. In our present study, the adopted constant limit of 0.50 ph cm -2 s _1 is larger than 
all individual limits in the catalog. 

All well defined sa mples above a given flux limit are subject to the Eddington effect 
(jEddingtonl Il913l . Il940l ). Random errors in the fluxes of individual objects if positive may 
cause them to become part of the sample, or, if negative, to be lost. Since there are generally 
more objects in the flux distribution below any given flux than above, the net effect will be 
that positive errors dominate in the sample. The effect on < V/Vmax > will be that for 
objects with a uniform distribution, the observed value will be larger than 0.5. 

We can evaluate the Eddington effect on < V/V max > of our sample through simulations. 
On the average, o for GUSBAD fluxes of the second brightest illuminated detector is 0.05 
ph cm -2 s _1 . Our adopted limit of 0.5 ph cm~ 2 s _1 (for detector 2) corresponds to ~ 1.21 
ph cm -2 s _1 for detectors 1 + 2, with o = 0.071 ph cm" 2 s" 1 . Simulations show that in this 
case the Eddington excess in < V/V mSQi > is 0.002, considerably less than the mean errors of 
< V/Vmax > for the spectral samples discussed below. 



2.2. Deriving Peak Energies 

Our goal is to derive the peak energy of the vF v spectrum. Given the low spectral 
resolution, we initially derive I(Ej) = EfN(Ei) for a representative energy E{ in each BATSE 
channel i, where N(Ei) is the photon flux density, in ph cm -2 s _1 keV -1 . We used the peak 
photon fluxes to derive the flux densities N(Ei) at E^ = 30, 70, 185 and 420 keV. For channels 
1-3, these energies are close to the geometric means of the energy band limits; for channel 
4, see below. 
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An initial estimate of the peak spectral energy is provided by the channel with the 
largest value of I(Ej). Table 1 shows for each of the four BATSE channels the number 
of GRBs having the peak at We can compare our -E p k,obs values with those given by 



Kaneko et al.l (120061 ) for peak fluxes. Their work was based on BATSE LAD data of 350 
bright GRBs in the BATSE catalogs. They used several LAD data types, and employed 
a number of spectral shapes beyond the Band spectrum. The Kaneko list has 219 GRBs 
that are in the GUSBAD catalog. The last column of Table 1 gives the average Kaneko 
peak energy for the GUSBAD sources in the Kaneko list. We adopt the Kaneko average 
for channel 4, since there is no well defined upper energy limit to its energy band. Only 
for channel 3 is a meaningful comparison possible: our result differs by only 10% from the 
Kaneko average. 

The averages of the Euclidean values of V/V mSuX for peak energies in channels 1-4 are 
0.44, 0.45, 0.32, and 0.33, respectively. It turns out that the transition from high to low 
values of < U/U max > takes place within channel 3. Since it contains a large number of peak 
energies, we subdivide those in channel 3 based on whether the ratio of I^E^) / I(E 2 ) is (a) 
less than 0.5, (b) in the range 0.5 — 1.0, or (c) larger than 1.0. Given that only 24 objects 
peak in channel 1, we combine channels 1 and 2. We end up with five spectral peak classes sp 
with < K/Knax > ranging from ~ 0.45 to ~ 0.30, see Table 2. The resulting -E p k,obs -V/14iax 
correlation is illustrated in Figure 1. 

Table 2 also gives the mean values of a 2 3, the photon spectral index derived from BATSE 
channels 2 an d 3, as given in the GUSBAD catalog. The a 23 — U/Vm ax correlation had been 



noted before (jSchmidtll200ll ) in a study similar to the present one. 



3. Derivation of Source Counts from the Luminosity Function 

Since we do not have redshifts for our sample of GUSBAD GRBs, we cannot derive the 
luminosity function directly. Instead we iterate the luminosity function by trial and error. In 
each step of the iteration we derive the predicted < V/V max > value for each of the spectral 
classes until they agree with those given in Table 1. 

We assume that the GRB luminosity function z, sp) of spectral class sp can be 
written as 

§(L,z,sp) = $ (L,sp)GR{z), (1) 

where L is the peak luminosity, &o(L,sp) is the z = luminosity function of class sp, and 
GR(z) the comoving GRB rate density, normalized at z — 0. We assume that <&o(L, sp) has 
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a gaussian distribution 

R (sp) -[log L -log L c (sp)] 2 
$o(£, sp) = 7 =— exp — 2 (2) 

with dispersion <T\ og L around a central peak luminosity L c (sp). Ro(sp) is the GRB rate at 
z = 0. 

Our approach will be to iterate the value of L c (sp) until the corresponding value of 
< V/V max >(sp) for the GUSBAD sample agrees with the observed value. This requires a 
full derivation of the expected source counts from the luminosity function. The role of &\ og L 
is primarily to make each spectral luminosity function sufficiently wide so that the total 
luminosity function is reasonably smooth. We ended up using <J\ og L = 0.5. 

In deriving fluxes from luminosities and redshifts, we employ the BATSE energy range 
(E 1 ,E 2 ) with Ei = 50 keV and E 2 = 300 keV. For an object at redshift z, the observed 
energy range (E 1 , E 2 ) originates in the range (E 1 (l + z), E 2 (l+z)) in the object's rest frame, 
whereas the luminosity refers to the range (Ei, E 2 ). The K-term is the ratio of the rest frame 
energies radiated in the two ranges, 

Klz) = . (3) 

J^EN(E,E pk ,a,(3)dE 

The Band photon spectrum is ususally described in terms of a break energy E . Here we 
use a Band spectr um N(E, E pk , a, (3) where E pk (sp) = (2 + a)E (sp) assuming that (3 < — 2 
(IBand et al.lll993l ). We adopt constant values of a = —0.8 and (3 = —2.6 for reasons that 
will be discussed in Section 5. 

The peak flux P(L, z) observed for a GRB of luminosity L at redshift z is 

P(L, z) = % K(z), (4) 

where (c/H )A(z) is the bolometric luminosity distance. We use the cosmological parameters 
H = 70 km s- 1 Mpc" 1 , Q M = 0.3, and Q A = 0.7. 

The integral peak flux distribution for GRBs of spectral class sp is, 

/pz(L,P,sp) 
$ (L,sp)dL J GR(z)(l + z)- 1 (dV(z)/dz) dz, (5) 

where z(L,P,sp) is derived from equation (4), V(z) is the comoving volume and the term 
(1 + z)~ l represents the time dilation. With this formulation, it is straightforward to derive 
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the differential source counts dN(> P,sp)/dP, as well as the average values of V/V ma , x , 
-Epk.obs, «23) etc. 

For a given rate function GR(z) the procedure to iterate the luminoxity function is as 
follows. Assume starting values for the central luminosity L c (sp) and the rest frame Band 
peak energy E p k(sp). The differential source counts together with V/V max = (P/ Pii m )~ 3 ' 2 
produce the expected values of < V/V max >(sp) for the GUSBAD catalog. Similarly, the 
expected values of -E p k,obs( s P) are obtained by weighting E pk (sp) / (1 + z) with the differential 
source counts. The iteration is repeated until the expected values of < V/V max >(sp) and 
-E P k,obs(sp) match the observed ones given in Table 2. 

It is worth noting that given GR(z) and the shape of the five spectral luminosity func- 
tions, the procedure leads for each sp to a single value for L c (sp) and E p k(sp). The primary 
unknown is the density function GR(z). We will use Swift data to test various forms of 
GR(z), see Sec. 5. 



4. Swift Data 



The Swift mission (IGehrels et al.ll2004l ) is in the process of detecting hundreds of GRBs. 
Thanks to an emphasis on rapid identification and communication, relatively many of these 
GRBs have observed redshifts. We use this data base to check predicted distributions of 
redshifts and luminosities based on luminosity functions with different GR(z). 

We have used Swift GRB data for the period Jan 1, 2005 - Sep 30, 2008. The observa- 
tions cover a field of view of 1.4 sr in the energy band 15 — 150 keV. We use observed peak 
fluxes over a one second time interval obtained from the High Energy Astrophysics Science 
Archive Research Center (HEASARC) and redshift information from J. Greiner's list 0. 
We derive isotropic-equivalent luminosities for the energy band 15 — 150 keV using the rest 
frame E pk — L iso correlation given in Section 6. In deriving these luminosities, we assumed 
a Malmquist correction AlogL = +0.36, see Section 6. 

The Swift m i ssion does not provide a flux limit above which the sample of GRBs is 
complete. iBandl (120061 ) has analyzed the sensitivity of Swift in detail. He concludes that 
the complexity of the trigger system maximizes the sensitivity, but "makes an accurate 
determination of this sensitivity at a given time very difficult if not impossible." Given these 
circumstances, we have to carry out an a posteriori estimate of the effective flux limit above 
which the burst list is complete. 



2 Available at http://www.mpe.mpg.de/~jcg/grbgen.html 
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Following the procedure outlined in Section 3, we can derive the predicted distributions 
of fluxes, luminosities and redshifts for Swift. The predicted flux distribution turns out to 
be virtually the same for all forms of GR(z) discussed in the next section. We show the 
distribution in Figure 2, together with the observed number of GRBs and the number of 
GRBs with observed redshifts. At large fluxes, the prediction appears to be around 9% 
below the observed numbers. The observed numbers start to deviate systematically below 
the prediction at P = 1.0 ph cm -2 s" 1 . We adopt this flux as the effective completeness 
limit for the Swift bursts. The total number of GRBs above the flux limit is 217, of which 
84 have redshifts. The fraction of GRBs with redshifts declines from ~ 0.6 at large flux to 
~ 0.35 at the limit. We will assume that a redshift fraction of 0.4 applies uniformly above 
the flux limit. 

With the Swift flux limit and the fraction with redshifts set, we are now in a position 
to derive from the luminosity function for a given GR(z) the predicted distributions of 
luminosity and redshift for the Swift sample. 



5. Exploring Luminosity Functions with Different GR(z) Functions 

As discussed at the end of Sec. 3, the only unknown in deriving the GRB luminosity 
function, given the -E p k,obs _ V/V^^correlation, is the redshift dependence GR(z). Therefore, 
we can view the excercise as one in which we explore which shape of GR(z) is compatible 
with the observations based on Swift data. 

We first explore the case where GR(z) is the galaxy star formation history (SFH). In the 
many published studies of star formation in the literature, the SFH rises fast with redshift 
by about an order of magnitude out to z ~ 1.5, beyond which it levels off; it tends to decline 
beyond z ~ 3. We use the analytical expression, normalized to z — 0, 

Psfh = (1.0 + (0.10/0.015)^/(1.0 + (^/3.4) 5 - 5 ), (6) 



for the SFH de scribed by lHopkins fc Beacoml (120061 ). based on an extensive compilation by 



Hopkinsl (120041 ) 



The predicted distributions of luminosity and redshift are compared to the Swift obser- 
vations in Figures 3 and 4. The predicted number of GRBs (with redshifts) is 79.6, while the 
observed number is 84. The prediction is ~ 5% below the observed number. It is obvious 
that the agreement with the observations is poor, both for the luminosities and the redshifts. 
We have carried out Monte Carlo simulations to find the probability P(logL > 51.5) that 
the observed number of luminosities above logL = 51.5 can be produced from the predicted 
distribution by chance. We find P(logL > 51.5) < 10~ 6 . Similarly, for redshifts above 3.0 
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the probability P(z > 3.0) < 10 6 , see Table 3. Even if we continue the SFH beyond the 
peak redshift z = 2.55 at its peak value without any downturn, the agreement remains poor, 
with P(logL > 51.5) < 10~ 6 and P(z > 3.0) = 0.007. Clearly, the SFH cannot represent 
GR(z). 

With little guidance as to what the shape of GR(z) could be, we use a simple schematic 
involving a (1 + z) power law rise, a plateau, and an exponential decline with z, as follows: 

GR(z) = (1 + z) m 0<z<z c (7) 

GR(z) = (1 + z c ) m z c <z<z d (8) 

GR(z) = (1 + z c ) m 10 k{z - Zd) z> z d (9) 

With the goal of producing redshift and luminosity distributions more in accord with the 
observed Swift distributions, we explored a number of combinations of the free parameters 
m and z c . We concluded from these explorations that to first order the ratio R = GR(z = 
4)/GR(z = 1.5) is crucial. This appears related to the fact that most Swift bursts have 
redshifts below 2. For the SFH considered above, the ratio R = 0.74. 

We present detailed results for two shapes of GR(z) that have R-values of 3.2 (A) and 
4.0 (B), respectively. Table 4 and Figure 5 illustrate the two cases. Case A rises by a factor 
of 32 to z = 3.0 and then remains constant. Case B rises slower to a plateau of 25 starting 
at z — 4.0. We have not actually included the exponential decline beyond z&. With ony two 
redshifts with z > 5.0 in the Swift sample, we have no observational leverage on the value of 
z d - 

The properties of the two models are given in Tables 3 and 4. Both are successful in 
producing luminosity and redshift distributions compatible with the Swift data, see Figures 
6 and 7. The redshift distributions predict that 3 — 4% of GRBs in the Swift sample have 
z > 6 and ~ 1% z > 8. Since we have allowed the density plateau to continue beyond 
z — 10, these percentages are overestimates. Our predictions are lower than most discussed 
in the literature. 

In deriving the luminosity functions, we found that the spectral indices a 23 given in 
Table 2 are best represented by the Band spectral parameters a = —0.8 and (3 = —2.6. We 
show the observed average spectral index a 2 3 for the five spectral classes versus -E p k,obs in 
Figure 8, together with the relation derived from the luminosity functions. 

We illustrate in Figure 9 the derivation of L c (sp) from the observed < V/V max > values 
for case A. The differences in the curves reflect the different values of E pk for the spectral 
classes. For the GUSBAD sample, the curves are only meaningful near the actual L c (sp) 
values. For, say, a deeper sample, the applicable parts would be at correspondingly higher 
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luminosities. 

Figure 10 shows the predicted distribution of redshifts for the GUSBAD sample, again 
based on case A. For sp — 1 essentially all the redshifts are expected to be below I. This 
component is very uncertain, as discussed below. The peak redshifts for the different spec- 
tral classes reflect their < V/V milx > values. This illustrates clearly that the Euclidean 
< V/V max > is a cosmological distance indicator. 

In Figure 11 we show the luminosity function for case A, as well as the luminosity 
distribution for the 1319 GUSBAD sources with P > 0.5 ph cur 2 s _1 . Also shown are the 
individual luminosity functions for the five spectral classes. The first peak of the luminosity 
function is contributed by spectral class 1. The lower half of its gaussian clearly plays no 
role, as it produces no objects in the luminosity distribution. The luminosity assigned to 
this class is uncertain, since the slope of the curve in Figure 9 is relatively shallow. Actually, 
if the < V/V max > for sp = 1 were only 1.2a larger, it could not be reproduced by any value 
of L c . Altogether, this suggests that the (large) z = density rates Ro for sp = 1 given in 
Table 4 are very uncertain. The second peak in the luminosity function is contributed by 
the large number of GRBs in spectral classes 2-5. Their combined z = rate is 0.09 — 0.22 
Gpc -3 y _1 , for models A and B, respectively. 



6. The -E p k,obs — L iso and E pk — L iso correlations 

.Epk correlations with radiated energy or luminosity are of interest in exploring the 
mechanism for the prompt emission of GRBs. They are also of practical interest in allowing 
an estimation of the redshift of GRBs with measured -E p k,obs- In this section, we discuss the 
derivation of the relevant isotropic-equivalent luminosity L iso , present the -E p k,obs — ^iso and 
-Epk — L iso correlations and briefly mention the problem of extracting individual redshifts 
from E pk correlations. 

As discussed in Section 3, the derivation of the luminosity function involves an iteration 
of the central luminosity L c (sp), where each spectral component is a gaussian with dispersion 
o"io g L- We chose a dispersion a\ og L= 0.5 that produces a reasonably smooth overall luminosity 
function. We want to use the L c (sp) luminosities in deriving the isotropic-equivalent peak 
luminosities L iso used in the correlations. 

It turns out, however, that L c (sp) varies considerably with u\ oz l. We explored using 
ciogL values ranging from 0.1 to 0.6 in deriving the luminosity function for cases A and B. 
For all values of a\ og L the agreement with the observed luminosities and redshifts of the Swift 
GRBs is good. The variation of \ogL c (sp) appears to be proportional to a\ ' L . 
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The explanation of the situation is as follows. As 0"i O gL is increased, the wings of 
the gaussian will provide more objects with higher and lower luminosities. The higher 
luminosities will be at larger distances and contribute lower V/V max values. The lower 
luminosities will contribute higher V/V max but the shift will be less due to the curvature of 
the relation between L c (sp) and V/V max , see Figure 9. The asymmetry results in a higher 
value of < V/Vjbbx > requiring a shift to a lower L c (sp) in order to fit the observed input value 
of < V/V max >. The reason why agreement with the Swift observations is not affected, is 
that the additional objects of lower luminosity contribute little to the observed distributions 
because they are observed over much smaller volumes. 

Given our understanding of the variation of L c (sp) with <J\ og L, it is clear that we want 
to use the unbiased values of log L c (sp) at <7i og L= in deriving the L iso . They are only 
0.01 — 0.03 larger than the logL c (sp) values for a\ og L= 0.1. The L iso values given in Table 4 
are 2 — 5 times larger than L c (sp). 

We now plot the L iso versus -E p k,obs to produce the ^pk.obs — -^iso correlation in Figure 12. 
The error bars for L iso reflect the effect of those in < V/V max >. The line drawn represents 
a regression of L iso on -E p k,obs, assuming that the -E p k,obs values are errorless, 



log Lu 



51.0 + 2.52l°f (log E } 



pk,obs 



2.08j£») 



(10) 



The derivation of the luminosity function involved an iteration not only of the central 
luminosity L c (sp), but also of the rest frame value of the Band peak energy E p ^(sp). Figure 
13 shows the E p k — L iso correlation. The line drawn represents a regression of L iso on E p k, 



log L h 



51.0 + 1.75l°;l(log£ pk - 2.36+°^ 



-0.42y 



The extraction of individual redshifts from the £pk,obs — Li SO correlation for sources with 
known £7 p k )0 bs is straightforward. Using the discrete E pkj0bs values of Table 2, we can produce 
redshifts for all 1319 GRBs without problems. Since the correlation uses the -E p k,obs i n the 
observer's frame, it can only be used for GRBs with a 50 — 300 keV flux larger than 0.5 ph 

— 2 —1 

cm s . 

The situation is entirely different for the E p ^ — L iso correlation. When we use the 
correlation to derive redshifts for sources with known -E p k,obs, we find that only 772 out of 
1319 GRBs, or 59% , yield a redshift. The situation is str i kingly simi lar to that repo r ted fo r 
the Amati relation jAmati et aDl2002[ ) bv lNakar k PirarJ pjoj and band k Preece! j2005h . 
who found that for a large fraction of BATSE bursts there is no redshift that satifies the 
relation. 



There are relativ ely few E^ — L; so correlations reported in the litterature (ISchaeferll2007l ; 
Yonetoku et al.ll2004l ). We will compare our correlation with the one given by Schaefer, based 
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on 62 GRBs with -E p k,obs values and redshifts, 

logL(bol) iso = 52.21 + 1.68(logE pk /300keV). (12) 

Here logL(bol) covers the energy band 1 — 10000 keV. Transformation of our 50 — 300 keV 
L iso values listed in Table 4 to 1 — 10000 keV increases them on the average by 0.61 in logL. 
Comparison of the bolometric luminosities yields 

< logL(Schaefer) - logL(present) >= +0.52 ± 0.22. (13) 



Next we consider the effect of Malmquist bias (IMalmquistll 1 936f ) , the difference in average 
luminosity of objects observed in a sample above a given flux limit and the average luminosity 
in space, 

A logL =< logL > obs - < logL > space . (14) 

In Euclidean space, the bias is proportional to cr 1 2 ogL . For cosmological objects with varying 
co-moving density, the situation is more complex. 

All studies of individual GRBs with redshifts aimed at deriving E pk correlations are 
subject to Malmquist bias since they can only be carried out above some flux limit. The 
Malmquist bias depends on the luminosity function (in space) of GRBs at given E pk . This 
function is not known at present. 

We do have clear evidence of Malmquist bias in the present study. Figure 11 shows that 
compared to the luminosity functions for the five spectral classes, the luminosity distributions 
predicted for the GUSBAD sample are shifted to considerably higher luminosities. The 
Malmquist biases for sp = 1 — 5 are 0.70, 0.53, 0.31, 0.28 and 0.36, respectively, for an average 
of 0.44. This value is tantalizingly close to the offset shown in eq. (13). However, our values 
are based on an assumed gaussian shape of the spectral luminosity functions with a dispersion 
of CT\o g L — 0.5, designed to make the overall luminosity function reasonably smooth. As long 
as we do not have reliable information about the distribution of log L iso at given E pk in space, 
we will not be able to assess the effect of Malmquist bias accurately. 



7. Discussion 

We present in this final section some commentary on the approach used in this paper 
and on the results. 

1. This work validates the use of < V/V max > as a cosmological distance indicator. It 
requires a complete sample of GRBs above a well defined flux limit. The GUSBAD catalog 
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with its uniform treatment of all sources is eminently useful for this purpose. The input in 
the form of < V/V max > forced us to a deductive approach where we used a given shape of the 
(spectral) luminosity function and iterated its center luminosity and rest frame E p k until the 
observed < V/V max > and -E p k,obs were fit. An important advantage of the deductive method 
is that the derivation of source counts from a given luminosity function can be carried out 
accurately. 

2. This approach is in contrast to the inductive method generally used when redshifts are 
available. This involves deriving densities in bins of redshift and luminosity and correlating 
these with redshift and/or luminosity to derive the luminosity function. In such an approach, 
it is essential to check that the resulting luminosity function produces the input source counts 
and < V/V max > correctly. If this is not t he case, predic tions about expected numbers of 
sources at large redshift will be unrealistic (jSchmidtl 2004 ). 



3. If GRBs were not subject to cosmological evolution, the input in terms of < V/V max >, 
-Spk.obs, and < «23 > would suffice to derive the luminosity function and hence predict source 
counts and the redshift and luminosity distributions. Comparing these predictions with the 
observations, such as from Swift, would then constitute a check on the cosmological model. 
In reality, there is evolution, so in the context of the present study the Swift observations of 
GRBs with redshifts essentially provide information about GR(z). 

4. We have indicated that our exploration of various forms for GR(z) indicate that the 
Swift observations seemed to require that the GR(z = 4 )/GR(z = 1.5) ratio be at least 3. 



None of the st ar formation r esults shown in Figure 1 of iHopkins fc Beacoml (120061 ) exhibit 
such a ratio. iLe fc Dermerl (120061 ) explored a variety of GRB rate curves and found the 
best fit for their model SFR6, for which GR(z = 4)/GR(z = 1.5) ~ 2. Given that GRBs 



tend to occur in blue, underluminous galaxies of limited chemical evolution (IFruchter et al. 



20061 ). it is not surprising that GR(z) is different from the main cosmic star formation, and 
concentrated to earlier cosmic times. A study of the hydrogen ionization rate, based on 



the Lyman-a forest (jFaucherl 120081 ) . concludes that the star formation rate is increasing for 
2 = 2 — 4, qualitatively much like our model B. 

5. We have been conservative in setting the effective completeness limit for the Swift 
data at P = 1.0 ph cm -2 s" 1 and using only data above the limit. It would be helpful 
if catalogs of Swift sources could include an indication for each source whether it can be 
considered to be part of a sample that is complete above a given photon limit over a given 
area. 

6. Even while not all individual GRBs can be fitted with a Band spectrum, it appears 
that for statistical work, the E p k — L iso correlation together with fixed values for a = —0.8 
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and (3 = —2.6 provide a surprisingly simple prescription. 

7. In comparing E p ^ — L lso and E p ^ — E iso correlations, the transformat ion from the 



actua l energy bands used, like 50 — 300 or 15 — 150 keV, to 1 — 10000 keV (lAmati et al. 



20021 ) introduces needless uncertainty. It would be preferable to report results in the observed 



energy bands. 

This research made use of data obtained from HEASARC, provided by NASA's Goddard 
Space Flight Center. It is a pleasure to thank Y. Kaneko for detailed information about 
spectra. 
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Table 1. Spectral Peak Energies 



Channel 


n a 


£ P k,obs(keV) a 


n b 


£ P k,obs(keV) b 


1 


24 


30 






2 


204 


70 


4 


75 


3 


608 


185 


99 


202 


4 


483 


420 c 


116 


417 



a GUSBAD sample 

feaneko et aD kood \ data 

c Value adopted from lKaneko et al 



120061) 
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Table 2. Peak Spectral Energy, < V/V max > and Photon Index 



Spec. CI. 


Channel 


n 


£ p k,obs(keV) 


< v/v max > 


< «23 > 


1 


1,2 


228 


65 


0.449±0.019 


-2.52 


2 


3 a 


185 


120 


0.395±0.021 


-1.77 


3 


3 b 


207 


175 


0.296±0.018 


-1.75 


4 


3 C 


216 


250 


0.283±0.018 


-1.44 


5 


4 


483 


420 


0.327±0.013 


-1.25 



a I(E 4 )/I(E 2 ) < 0.5 

b 0.5 < I(E 4 )/I(E 2 ) < 1.0 

c 1.0 < I{E A )/I(E 2 ) 
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Table 3. Comparison of Luminosity Function Models 



Model 


SFH 


A 


B 


GR{z) 


PSFH a 


(1 + z) 2 - 5 


(1 + z) 2 - 


z c 




3.0 


4.0 


Zd 




10.0 


10.0 


k 




0.0 


0.0 


P(logL > 51.5) b 


< io~ 6 


0.16 


0.18 


P(z > 3.0) c 


< KT 6 


0.93 


0.97 


f(z > 2) d 


0.198 


0.448 


0.426 


f(z > 4) d 


0.009 


0.126 


0.159 


f{z > 6) d 


0.000 


0.037 


0.049 


f(z > 8) d 


0.000 


0.012 


0.017 



a Analyti cal expression for the sta r formation his- 
tory SFH (IHopkins fc Beacomll2006l ). see eq. (6) 

Probability that the model produces the ob- 
served number of Swift bursts with logL > 51.5 

"^Probability that the model produces the observed 
number of Swift bursts with z > 3 

d Fraction of Swift bursts with P > 1.0 ph cm -2 
s -1 above the given redshift 
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Table 4. Spectral Class Luminosity Functions 







Case A 






Case B 




Spec. CI. 


logL c a 


logL iso b E C 


Ro d 


logL c a 


logL iso b E C 


R d 



1 


49.1118$ 


49.8ll°| 


73 


6.8 C 


48.9318$ 


49.63181? 


69 


13.0 C 


2 


50.84±8;|| 


5 1 . 53 q 


240 


0.048 


50.5418-22 


5i.i4i8:21 


203 


0.149 


3 


51.67jlg;il 


51.97t°£ 


506 


0.0089 


5i.7H8ii 


52.0718:8? 


518 


0.0133 


4 


51.70lg;S 


ci nc+0.06 

oi.yo_ 06 


762 


0.0078 


ci 7 c+0.12 
oi - ' °-0.13 


c 2 nc+0.06 

oz.uo_ 006 


792 


0.0112 


5 


CI oq+0.08 


51 7fi+ - 04 

Ji - ' u -0.04 


1217 


0.0225 


ci /fo+0.08 


ci oc+0.04 

51.85_ 004 


1244 


0.0336 



a Central isotropic-equivalent luminosity of the gaussian in the 50 — 300 keV energy band 
for a\ ogL = 0.5, in erg s" 1 . Errors correspond to those for < V/V max >, see Table 2. 

b Isotropic-equivalent peak luminosity in the 50 — 300 keV energy band, in erg s -1 , see Sec. 
6. Errors correspond to those for < V/V max >, see Table 2. 

c Band spectrum break energy in the rest frame, in keV 

d GRB density rate at z = 0, in Gpc~ 3 y _1 

c Rates for sp — 1 are very uncertain, see Sec. 5. 
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100 



200 



500 



E 



pk.obs 



(keV) 



Fig. I. — < V/V max > values versus the observed peak spectral energy E^obs based on 1319 
GUSBAD GRBs with P > 0.5 ph cm" 



2 s" 1 . 
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SWIFT Differential Source Counts 

o 

o 




P (ph cm 2 s ') 

Fig. 2. — Differential source counts for Swift, observed number (dots), those with redshifts 
(circles) and the predicted number based on a luminosity function (line). Based on this plot, 
we estimate that the effective completeness limit for Swift bursts is P = 1.0 ph cur 2 s -1 . 
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Fig. 3. — Luminosity distribution of Swift sources: observed (dots) and predicted (curve) if 
GR(z) equals SFH, the typical galaxy star formation rate, see .eq. 6. 
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Fig. 4. — Redshift distribution of Swift sources: observed (dots) and predicted (curve) if 
GR(z) equals SFH, the typical galaxy star formation rate, see .eq. 6. 
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Redshift z 

Fig. 5. — Normalized GRB rates GR(z). The typical star formation rate SFH produces too 
few redshifts larger than 2. Models A and B predict luminosity and redshift distributions 
that are compatible with the Swift observations. 
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L peak (l5_15 ° keV) (erg SeC_1) 

Fig. 6. — Luminosity distribution of Swift sources: observed (dots) and predicted for models 
A (full line) and B (dashed line). 
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Redshift z 



Fig. 7. — Redshift distribution of Swift sources: observed (dots) and predicted for models A 
(full line) and B (dashed line). 
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100 200 500 

E pk> ob S (keV) 

Fig. 8. — Observed value of the photon index «23 versus the peak spectral energy -E p k,obs 
(points) and average values derived for models A and B (curve) with a = —0.8 and (3 = —2.6. 
The error bars for most of the points are too small to show. 
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T 1 1 1 1 — I I I | 1 1 1 1 1 — I I I | 1 1 1 1 1 — I - TT 




I I I I I l_l I I I I I I I l_l I I I I I I I l_l I I I I I ■ I ■ 

10 49 10 50 10 51 10 ! 

L c (50-300 keV) (erg s" 1 ) 

Fig. 9. — Illustrating the derivation of L c (sp) from the observed values of < V/V max >. The 
points correspond to model A, see Table 4. 
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Redshift z 

Fig. 10. — The redshift distribution for the GUSBAD sources, for each of the five spectral 
classes, based on the luminosity function for model A. 
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Fig. 11. — Left side: the luminosity function for model A (full line) and the five spectral class 
luminosity functions (dashed lines). Right side: the corresponding luminosity distribution 
ra(logL) in steps of A log L = 0.1 for 1319 GUSBAD sources with P > 0.5 ph cm" 2 s -1 . 
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Fig. 12. — The isotropic-equivalent luminosity L iso versus the observed -E p k,obs for the five 
spectral classes, for models A (circles) and B (triangles). 
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Fig. 13. — The isotropic-equivalent luminosity L iso versus the rest frame E p ^ for the five 
spectral classes, for models A (circles) and B (triangles). 



